Quantization of a Free Particle Interacting Linearly with a Harmonic Oscillator 
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We investigate the quantization of a free particle coupled linearly to a harmonic oscillator. This 
system, whose classical counterpart has clearly separated regular and chaotic regions, provides an 
ideal framework for studying the quantization of mixed systems. We identify key signatures of the 
classically chaotic and regular portions in the quantum system by constructing Husimi distributions 
and investigating avoided level crossings of eigenvalues as functions of the strength and range of 
the interaction between the system's two components. We show, in particular, that the Husimi 
structure becomes mixed and delocalized as the classical dynamics becomes more chaotic. 
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Typical classical Hamiltonians systems are nei- 
ther fully integrable nor fully chaotic but instead 
possess mixed dynamics, with islands of stabil- 
ity situated in a chaotic sea. In this paper, we 
investigate the quantization of a recently-studied 
system with mixed dynamics |l2l |. This example 
consists of a free particle that moves around a ring 
that is divided into two regions. At the bound- 
aries between these regions, the particle is kicked 
impulsively by a harmonic oscillator (in a manner 
that conserves the system's total energy), but the 
particle and oscillator otherwise evolve freely. Al- 
though the system is not generic, its separation 
into regular and chaotic components also allows 
more precise investigations (both classically and 
quantum-mechanically) than is typically possible, 
making this an ideal example to achieve a bet- 
ter understanding of the quantization of mixed 
systems. By examining avoided level crossings 
and Husimi distributions in the quantum system, 
we investigate the quantum signatures of mixed 
dynamics, demonstrating that the Husimi struc- 
tures of nearby states become mixed and delocal- 
ized as chaos becomes a more prominent feature 
in the classical phase space. 



INTRODUCTION 

Investigations of the quantization of chaotic systems 
have become increasingly prevalent as physicists conduct 
more experiments at small scales and design an increas- 
ing number of devices that exploit the physics at such 
scales [HI, 0, EH EH Experiments on quantum 

chaos, conducted using microwave cavities [3 El], atom 



optics [HEH, and other systems, have examined phenom- 
ena that are both fundamental and diverse-ranging from 
the decay of quantum correlations [l[ to localization in 



quantum wave functions [28| and chaotic scattering 31 1 . 
Despite this wealth of research, however, it is still not 
entirely clear how to understand the notion of chaos in 
quantum mechanics. Quantum wave functions satisfy a 
linear differential equation (the Schrodinger equation), 
so sensitive dependence on initial conditions and the ex- 
ponential divergence of nearby trajectories-key compo- 
nents for defining classical chaos-cannot be used to define 
quantum chaos. Nevertheless, quantum analogs of clas- 
sically chaotic systems do possess identifying features, so 
the quantizations of chaotic systems can be distinguished 
from the quantizations of integrable (regular) ones. 

Typical classical Hamiltonian systems are neither fully 
chaotic nor fully regular; rather, they have "mixed" dy- 
namics (i.e., a divided phase space), with islands of sta- 
bility ("KAM islands") situated in a chaotic sea. Be- 
cause generic mixed systems are very difficult to ana- 
lyze, there have been numerous attempts to construct 
Hamiltonian systems with mixed dynamics that allow an 
exact, rigorous analysis. Previously studied examples in- 
clude billiards [1, E|, M. l38l Ep- Fermi accelerators and 
bouncin g-ba ll models [13l. 29 , 43| , and kicked rotors and 

tops [laMEEE!- 

In the present paper, we investigate the quantization of 
a one-dimensional free particle interacting linearly with 
a one-dimensional harmonic oscillator. This recently- 
studied example has mixed dynamics with well-separated 
integrable and chaotic regions The clean separa- 

tion between different types of behavior helps simplify 
comparisons between the dynamics of the classical sys- 
tem and that of its quantization and makes this system 
a very illuminating one for studying the quantization 
of systems with mixed dynamics. The investigation of 
classical-quantum correspondences is extremely difficult 
for generic mixed systems, which possess an infinite hier- 
archy of KAM islands and intricately mixed chaotic and 
integrable regions. This makes the study of identifying 
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FIG. 1: [Color online] The configuration space of the particle. 
Its position on the ring is denoted q. 



features of chaos in carefully-chosen examples particu- 
larly important. 

The rest of this paper is organized as follows. First, 
we briefly review the classical system studied in [l2j ]. 
We then quantize this system and examine its sharp 
and broad avoided level crossings as the relative length 
of the interaction versus non-interaction region is var- 
ied. We illustrate our observations using Husimi distri- 
butions, which also allow a comparison with the classical 
dynamics. Finally, we summarize our results and present 
additional technical details of our investigation in two 
appendices. 



THE CLASSICAL SYSTEM 



Motivated by investigations of electron-phonon inter- 
actions m_ condensed matter physics d, 0, [H, [H, [H, 



261 . 1271 1341 l45j . De Bievre, Parris, and Silvius recently 



performed an analysis of a closed (classical) Hamilto- 
nian system consisting of an interacting one-dimensional 
free-particle and a one-dimensional harmonic oscillator 
[l2T ] . The example they investigated consists of a classical 
particle moving on a ring divided into two sections (see 
Fig. [I]): a region of length 2 called the "interaction re- 
gion" and one of length L called the "non-interaction" (or 
"uncoupled") region. At the boundaries between these 
regions, the particle is kicked impulsively by a harmonic 
oscillator (which moves on a line rather than on the ring) , 
but the particle and oscillator are otherwise uncoupled. 



However, the harmonic oscillator oscillates about differ- 
ent equilibrium points in the two regions; as discussed 
below, the equilibrium position in the interaction region 
depends on the strength of the interaction. 

The (non-dimcnsionalizcd) Hamiltonian describing 
this system is 
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(1) 



where p and II are, respectively, the particle and oscilla- 
tor momenta, q and <!> are the particle and oscillator po- 
sitions, a describes the strength of the particle-oscillator 
interaction, and x(l) takes the value 1 in the interaction 
region and everwhere else. We choose coordinates so 
that the interaction and non-interaction regions occur, 
respectively, when q g [0, 2] and q £ (2, 2 + L), where L 
is the length of the non-interaction region and 2 + L is 
identified with 0. The only system parameters that can 
be varied are a and L. 

Let's review some of the basic qualitative dynamics 
of JT]) 12]. The system achieves its ground state when 
the particle is in the interaction region and the oscilla- 
tor and particle are both at rest. When transitioning 
between the interaction and non-interaction regions, the 
particle's momentum changes discontinuously and its po- 
sition changes continuously. When the particle enters the 
interaction region, the harmonic oscillator experiences an 
interaction force that shifts its equilibrium position from 
to a. (The oscillator's momentum and amplitude de- 
pend continuously on time.) One uses conservation of en- 
ergy to compute the impulse that the oscillator imparts 
to the particle at the transition points. See Ref. [12j for 
further details. 

The uncoupled Hamiltonian, given by equation ([1]) 
with a = 0, is integrable. It possesses two independent 
integrals that correspond, respectively, to the energy of 
the harmonic oscillator and the momentum of the par- 
ticle. The latter integral leads to an SO(2) rotational 
symmetry in the particle's configuration space. Indeed, 
§Q is invariant under symplectic transformations 

T A : (p, q, n, $) i * (det(A)p, Aq, U, $), A e 0(2) , (2) 

forming a symmetry group isomorphic to 0(2). However, 
for a ^ (i.e., the generic case), the Hamiltonian (flj no 
longer possesses the two integrals and is only invariant 
under the subset of the symplectic transformations ([2]) 
with either A = 1 or A a reflection about the line pass- 
ing through q = 1 and q = 1 + L/2 (see Fig. [1]). Such 
transformations form a subgroup of 0(2) isomorphic to 
Z2 = {1,-1}, reducing the symmetry of the system to a 
parity symmetry about an axis of the ring in Fig. [T] 

A particularly interesting facet of this system is the 
clean separation of the integrable and chaotic regions in 
its phase space. Phase portraits of the system possess 
two characteristic integrable regions (among other struc- 
tures). The first, which exists for system energies varying 
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Degeneracy Lifting of the 17 lh -20 lh Levels for L=2 



First 100 Energy Levels for L=2 




FIG. 2: [Color online] Lifting of the double degeneracy in the 
eigenenergies for a > 0. 



from the ground-state energy to a critical positive energy, 
arises when the particle never leaves the interaction area. 
For small positive energies outside this integrable region, 
the motion appears to be fully chaotic without any ad- 
ditional KAM structures near the boundary between the 
two regions [l2l |. The second integrable region is an el- 
liptic KAM island centered on the equilibrium point that 
arises from the orbit in which the particle traverses each 
section of the ring exactly once per period. More gener- 
ally, the distinction between chaotic and integrable dy- 
namics in this system is clear in a variety of situations. 



THE QUANTUM SYSTEM 



We quantize Eq. JT]) using canonical quantization 11 1 . 
In so doing, we assume that the particle and harmonic os- 
cillator act as bosons with no internal degrees-of-freedom 
(DOF) and impose the following commutation relations: 



\p, q ] = -i, [p,n] = o, [n,$] = 
[«,*] = o, M] = o, [n,<?] = o. 

With the coordinate-space identifications 

.d .d 



we obtain the quantum Hamiltonian 

d 2 



(3) 



(4) 



(5) 



For the uncoupled system [Eq. ([5|) with a = 0], the 
time-independent Schrodinger equation H\tp) = E\ip) is 
separable, so one just needs to determine the eigenstates 
of the one-dimensional harmonic oscillator (moving on a 
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FIG. 3: [Color online] (a) The first 100 energy levels as a 
function of the interaction strength a for Eq. ([5]) with a non- 
interaction region of length L — 2. (b) Magnification of the 
avoided crossing between the 18th and 19th levels from panel 
(a). The inset shows a further magnification, and the la- 
bels designate where we calculated Husimi distributions (see 
Fig. EJ. 



line) and the free particle confined to a ring of length 
2 + L as separate problems (both of which admit closed- 
form solutions). Let {\ip^ art )} n and {\ip^. sc }}k denote 
eigenstates of the particle and the harmonic oscillator, 
respectively, so that {\^^ rt ) <8> |V , fe SC )}n,fc are eigenstates 
for the uncoupled system. In fact, these states form a 
basis for the Hilbert space of either the coupled or un- 
coupled system. We represent the Hamiltonian ([5|) as 
an infinite matrix using this basis (see Appendix I) and 
approximate its eigenvalues and eigenstates using those 
of a truncation of the matrix. Because the eigenenergies 
of the particle are doubly degenerate, the eigenenergies 
of the full uncoupled system are also doubly degenerate. 
Second-order perturbations in a lift this degeneracy for 
a > and L / 2; numerical calculations indicate that 
this degeneracy is also lifted for L = 2 (see Fig. [5]). 
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First 100 Energy Levels for L=12 
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Cluster of Avoided Crossings between the 119 -121 s Energy Levels for L=12 



11.52 



Energy Levels for L=12 




(b) 

FIG. 4: [Color online] (a) The first 100 energy levels as a func- 
tion of a for Eq. ([5]) with a non-interaction region of length 
L — 12. (b) Magnification of panel (a) illustrating broad 
avoided crossings. 



The uncoupled quantum Hamiltonian commutes with 
both the particle momentum p and the Hamiltonian 
H sc — (n 2 +<!> 2 )/2 describing an isolated harmonic oscil- 
lator; these are quantized versions of the two independent 
integrals of the classical Hamiltonian. The uncoupled 
quantum Hamiltonian is thus integrable [Blj . l59j . The 
symmetry-breaking of the classical Hamiltonian is mir- 
rored in the quantum system, as [H,p] — a$>[x(q),p] ^ 
and [H, H osc ] = ax(q)[U 2 , $]/2 ^ for a > and generic 
values of L. The parity symmetry in q and the time- 
reversal symmetry, described for the coupled classical 
Hamiltonian, are still present in (|5|), but there are no 
obvious continuous symmetries for a > 0. 



AVOIDED CROSSINGS 

As the coupling parameter a is varied, the eigenval- 
ues of ([5]) can approach each other very closely or even 
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FIG. 5: [Color online] An isolated cluster of two sharp avoided 
crossings and a broad avoided crossing between the 119th, 
120th, and 121st levels for L = 12. The incoming slope of 
the 119th level is imparted to the outgoing slope of the 121st. 
The labels designate locations at which we calculated Husimi 
distributions (see Fig. I10p . 



cross. If the Hamiltonian is invariant under a symmetry 
transformation for a certain range of a, it can be block- 
diagonalized by exploiting this symmetry. One does this 
by choosing each block to be invariant under the symme- 
try transformation. Energy levels belonging to different 
blocks can cross as a is varied 

US 53, Hi- ° 

n the other 

hand, if a quantum Hamiltonian has no symmetries other 
than time reversal, then such a level crossing is called an 
"accidental degeneracy" and requires the confluence of 
two parameters jZJ]. In this case, most levels that ap- 
proach each other end up avoiding one another instead 
of crossing. 

Classical chaotic systems have fewer constants of mo- 
tion than DOF and thus have fewer symmetries than in- 
tegrable systems with the same number of DOF. One ex- 
pects the quantization of these two situations to exhibit 
signatures of this difference 58, 59] . Hence, the quantiza- 
tion of a chaotic system should possess fewer level cross- 
ings than the quantization of an integrable one, and the 
presence of numerous avoided crossings between energy 
levels provide a signature of chaotic regions in the clas- 
sical system. Indeed, avoided level crossings ar e typica l 
features of quantum chaotic systems 17. 24. 4fj| . 51 , 5f| . 

Figure^ shows the first 100 energy levels as a func- 
tion of a in a system with L = 2 using an a-step size 
of 1.5 x 10~ 3 and a 2916 x 2916 truncated Hamiltonian 
matrix. (The a-step size is the distance between suc- 
cessive values of a for which we calculate eigenvalues 
and eigenvectors.) One observes a multitude of appar- 
ent level crossings. Refining the numerical computation 
at some of these apparent crossings shows that they are 
actually avoided crossings at which the slopes of the en- 
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ergy level curves are exchanged. Avoided crossing 



s of this 



nature are known as "sharp" avoided crossings [4l|, [51 1 . 
In passing through such avoided crossings, the participat- 
ing levels exchange their eigenstate structures, behaving 
as though they had entered a level crossing [53|. As we 
discuss in more detail later, we have verified numerically 
that this indeed occurs for our system. Similar phenom- 
ena have also been observed in other systems, such as a 
sinusoidally driven particle in a square potential well [5ll ] 
and a hydrogen atom in a strong magnetic field [13] • 

Figure [3Jd shows a magnification of an avoided crossing 
between the 18th and 19th levels from Fig. [3^ using a 
refined a-step size of 2 x 10~ 6 . In general, the maximum 
a-step size at which the avoided crossings in Fig. [3^ can 
be resolved is 0(1O -6 ). As a result, it is time-consuming 
to verify numerically that all of the apparent crossings 
are actually very sharp avoided crossings. Although we 
have only observed avoided crossings, we have not ruled 
out that actual level crossings might occur. Using parity 
symmetry, the Hamiltonian ([5]) can be separated into the 
direct sum of two blocks. As a is varied for a fixed value 
of L, energy levels from different blocks might cross. 

Figure UK shows the first 100 energy-level curves for 
L = 12. Our numerical computations, for which we used 
a 4096 x 4096 truncated Hamiltonian matrix, verify that 
the frequency of avoided crossings increases as L grows. 
Even more interesting are the broad avoided crossings, 
which become more prevalent as a increases. Figure [4Jd 
shows energy-level curves for a £ [5.8,7.2]. The a-step 
size for which the sharp avoided crossings in this figure 
can be resolved is O(10~ 3 ), about one thousand times 
larger than that required to resolve the sharp avoided 
crossings in Fig. [3^- We have also observed that the 
prevalence of broad avoided crossings increases as L in- 
creases. In general, broad avoided crossings tend to occur 
in nearly isolated clusters in which only a subset of the 
initial slopes of the participating levels are exchanged af- 
ter the sequence of crossings in the cluster |60| . This be- 
havior is caused by the presence of "overlapping" avoided 
crossings in the cluster (by "overlapping," we mean that 
there are ranges of a values for which multiple avoided 
crossings in the cluster are taking place). Such clusters 
usually induce nontrivial structural exchanges between 
the participating eigenstates and are suspected to be an 
indication of chaos 
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521 ]. Indeed, prior work on a si- 
nusoidally driven particle in a square potential well |5l| 
has indicated that broad avoided crossings produce su- 
perpositions of eigenstate structure rather than complete 
exchanges. We will use the term "mixing" to refer to such 
superpositions. 

An example of an isolated cluster of avoided cross- 
ings is shown in Fig. [5] This cluster occurs for L = 12 
and shows the 119th-121st levels as they experience two 
sharp avoided crossings and one broad one. A broad 
avoided crossing begins between the 120th and 121st lev- 
els, but as a is increased, two sharp avoided crossings 



(a) 
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FIG. 6: [Color online] (a) Husimi distribution localized 
around quasiperiodic orbits in a KAM island. Lighter regions 
have higher probabilities and black regions are ones with zero 
probability. The classical SOS is plotted in turquoise. To 
facilitate the comparison between the quantum and classi- 
cal dynamics, we include plots of a few of the orbits in the 
integrable region, (b) Husimi distribution localized around 
multiple KAM islands. 



(a) 




(b) 




FIG. 7: [Color online] (a) Husimi distribution delocalized 
throughout the chaotic sea on the left half of the plot, (b) 
Husimi distribution and classical SOS for the chaotic eigen- 
state in (a) taken on the surface defined byg = 2 + i/4 = 12 
(one fourth of the way into the non-interaction region). 



cause structural exchanges that induce a broad avoided 
crossing between the 119th and 120th levels. The figure 
gives the appearance of a broad avoided crossing over- 
lapping with the second sharp avoided crossing. Con- 
sequently, the slopes of the participating levels in the 
second sharp avoided crossing are not completely ex- 
changed. We have observed that such behavior often 
arises in avoided crossings (one of which is almost always 
broad) that occur simultaneously as a is varied. While 
the complex interactions between energy levels (as a func- 
tion of a) may make it impossible to precisely identify 
broad versus sharp avoided crossings in isolated clusters, 
approximate distinctions such as that just discussed al- 
low one to classify these interactions and determine their 
role in structural exchanges between eigenstates (see the 
discussion below). 

One can gain insight into the energy level dynamics 
for highly excited states (i.e., in the deep semiclassi- 
cal regime) using level spacing statistics O, IH| . Once 
desymmetrized, completely integrable systems obey Pois- 
son statistics and completely chaotic ones obey Wigner 
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statistics. The quantization of mixed systems with well- 
separated integrable and chaotic regions (such as the 
present one) have energy level statistics that interpo- 
late between these two extremes, with a distribution of a 
precise form that is conjectured to follow Berry-Robnik 
statistics [f| Ho|, l39| - In the present paper, we focus on 
eigenstate structures rather than semiclassical dynamics, 
and we accordingly investigate Husimi distributions be- 
low and leave the analysis of energy level statistics for a 
future publication. 



HUSIMI DISTRIBUTIONS 

Although there is no equivalent of classical phase space 
trajectories in quantum mechanics, there are suitable 
analogs. In particular, the Husimi distribution is often 
used in the study of quantum chaos ii. Given a 
quantum state \ip), its Husimi distribution H^{p,q) is 
defined by the projection of \ip) onto a coherent state 
\ip(p, q )) localized around (p,q): H^,{p,q) cx | \ip) \ 2 ■ 
For a system with Euclidean topology, a coherent state 
localized at (p, q) is a Gaussian with position-space rep- 
resentation localized around q and momentum-space rep- 
resentation localized around p. The system ([5]) possesses 
a cylindrical phase-space topology, as the particle posi- 
tion q is a periodic variable and the momentum pel. 
In Appendix II, we construct the coherent state for this 
topology from the Euclidean coherent state |if|. 

Coherent states provide excellent quantum analogs of 
classical particles when visualized as wave packets that 
minimize the position- momentum uncertainty product. 
The projection onto these particle-like states can thus be 
viewed intuitively as a sort of classical smearing. One 
then interprets the Husimi distribution as a probability 
distribution in phase space 0, allowing one to under- 
stand the dynamics of a quantum system in an analogous 
manner to p hase portraits (in particular, surfaces of sec- 
tion 29|, Iff) of its classical counterpart. See 0,0, 
for additional discussions and applications. For the clas- 
sical system ((TJ), we plot Poincare surfaces of section 
(SOS), for p > and various values of q (often q = 1), 
on the projection of phase space onto the harmonic os- 
cillator's coordinates [i.e., on ($, n)-space] [l2j. As dis- 
cussed in Appendix II, we calculate Husimi distributions 
for eigenstates of the Hamiltonian and project them 
onto the quantum oscillator phase space. Using the eigen- 
state energy and the same values for the system param- 
eters a and L, we compare the Husimi distributions to 
the corresponding classical SOS. 



(a) 



(b) 



Localization Around Classical Features 

To investigate the signatures of chaotic and integrable 
dynamics on the quantum system, we compare quan- 



FIG. 8: [Color online] (a) Husimi distribution and classical 
SOS for the chaotic eigenstate in Fig.[7£i taken on the surface 
defined by q = 2 + L/2 = 22 (in the middle of the non- 
interaction region), (b) Husimi distribution showing both 
chaotic and integrable features while having only one con- 
nected component. The structure is relatively delocalized 
throughout the enclosed area on the left side of the available 
phase space. 



turn Husimi structures to corresponding classical SOS 
0, E3, [S3, EH - As expected, eigenstates sometimes show 
localization around distinctive features of an SOS. A 
prominent manifestation that often occurs is a strong lo- 
calization around KAM islands, which is sometimes ac- 
companied with Husimi density throughout the chaotic 
sea. Figure [6] shows two such states: In panel (a), the 
Husimi structure is strongly localized around quasiperi- 
odic orbits in an integrable region of the classical SOS; 
in panel (b), the structure is strongly localized in and 
around multiple KAM islands. In this figure, we dis- 
play the Husimi distribution as a contour plot that varies 
from black (zero probability) to white (about 10 -3 ), with 
lighter regions having higher probability than darker re- 
gions. We overlay the classical SOS in turquoise. We 
observe that the quantum manifestation of a classical in- 
tegrable region is a Husimi structure strongly localized 
inside the region. By analogy with the classical dynam- 
ics, we refer to such states as "regular." 

Figure [7^, shows an eigenstate whose Husimi structure 
is located in the chaotic sea of its corresponding classi- 
cal SOS (particularly around the edge on the left half of 
the SOS). We observe derealization of the distribution 
throughout the chaotic sea, providing a signature of the 
predominantly chaotic dynamics of the classical SOS. By 
analogy with the classical dynamics, we refer to quan- 
tum states with such structures as "chaotic." For Husimi 
distributions taken at q = 1, we have only observed 
structures localized around the left half of the available 
phase space. However, distributions constructed for q 
values in the non-interaction region of phase space have 
structures with positive density throughout the avail- 
able phase space. Figures [TJd and[8ji, for example, show 
Husimi distributions (and their corresponding classical 
SOS) for the chaotic state in Fig. [7K taken on the sur- 
faces defined by q = 2 + L/4 = 12 (a fourth of the way 
into the non-interaction region) and q = 2 + L/2 = 22 
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(the middle of the non-interaction region), respectively. 
These Husimi distributions still display a strong localiza- 
tion on the left but nevertheless have a significant density 
throughout the available phase space. 

The enhanced density around the left edge of the 
chaotic sea is present for all Husimi distributions with 
chaotic structures. We observe this as well for Husimi 
distributions taken in the non-interaction region. As il- 
lustrated in Fig. [SJd, this occurs even for large energies, 
for which the chaotic region in the left half of the SOS 
retreats to the edge of the available phase space. The 
location of this enhanced density coincides with the left 
boundary of the intersection between the available os- 
cillator phase space in the interaction region and the 
available oscillator phase space in the non-interaction re- 
gion. As discussed in detail in Ref. [13], the portion 
of oscillator phase space corresponding to the interac- 
tion region is given by a disk-shaped region of radius 



V2E- 



centered at (en, 0) and that corresponding to 



the non-interaction region is given by a disk-shaped re- 
gion of radius \J2E centered at the origin. Comparing 
Figs. [7^ and[7J}, we observe that this intersection in clas- 
sical phase space occurs for $ w —2 between n w —2 and 
IT w 2. In both plots, the Husimi density is clearly en- 
hanced in this area, in correspondence with the dynamics 
of the classical system. 

In Fig.[SJ), we show a Husimi distribution that contains 
a mixture of regular and chaotic features. Although the 
regular and chaotic structures appear to be disconnected 
at first glance, one can see upon very close inspection that 
they are actually connected by small "bridges" with pos- 
itive Husimi density. Such bridges appear more promi- 
nently between the regular structures in Fig. [BJd. The 
bridges that are often observed between Husimi struc- 
tures localized in different classical phase space regions 
(for example, between integrable and chaotic regions) 
serve as channels through which density continuously 
flows as a is varied through an avoided crossing. We did 
not observe any Husimi distributions in which the regular 
and chaotic components were completely disconnected. 

In concluding this subsection, we remark that there is a 
considerable body of work on "regular" versus "chaotic" 
eigenfunctions. According to Percival's conjecture, eigen- 
functions localize either on integrable or chaotic regions 
of the underlying classical phase space, so that the par- 
tial level density of regular states is given approximately 
by the fraction of volume of the integrable region rel- 
ative to the available phase space [!, [j, [36j] . A recent 
paper by Marklof and O'Keefe [HI contains rigorous re- 
sults on such an extension of quantum ergodicity theory 
to a general class of quantum unitary maps whose un- 
derlying classical system has a divided phase space. It 
would be extremely interesting to obtain similar results 
for mixed systems like the present one. 



Exchange and Mixing of Husimi Structures at 
Avoided Crossings 

Sharp Avoided Crossings 

In Fig. [SI we depict Husimi distributions for the 18th 
and 19th energy levels of ([5]) with L = 2. This plot re- 
veals the structural changes that occur as the two levels 
encounter the avoided crossing in Fig. [3)3. In the top 
panels, one can see that the 18th level is a regular state, 
whereas the 19th level is chaotic. The middle panels 
give snapshots near the closest point of the encounter. 
Here, the Husimi distributions appear as mixtures of the 
two initial distributions. Additionally, the regular and 
chaotic portions of the structure are connected. As a is 
varied, Husimi probability flows continuously between in- 
tegrable and chaotic regions. The bottom panels, depict- 
ing snapshots from long after the encounter, show that 
the two levels have completely exchanged their struc- 
ture through the avoided crossing, leaving the aggregate 
Husimi structure unchanged. This provides an example 
of a smooth exchange of character in a sharp avoided 
crossing, which has also been observed in other quan- 
tum chaotic systems 2l|, [5lj . The avoided crossings that 
we have observed between chaotic and regular states and 
between two chaotic states have all been sharp ones. 



Broad Avoided Crossings 

As described previously, broad avoided crossings in the 
energy-level curves typically occur in nearly isolated clus- 
ters in which the initial slopes of the curves are not fully 
exchanged after the sequence of avoided crossings. This 
leads to a "mixing" of Husimi structures rather than a 
complete exchange [Hl| , which we have verified is the case 
for the Hamiltonian ([5|). Figure [TOl shows an example of 
such a mixing in the cluster of avoided crossings from the 
119th through the 121st energy levels from Fig. [5[ The 
localization of the initial Husimi structure of the 119th 
level indicates that this eigenstate is regular. The 120th 
level is localized near an integrable region, although there 
is a significant Husimi density in the chaotic sea, and 
the 121st level displays a strong chaotic localization with 
a slight density near the integrable region. The 119th 
and 120th levels leave the avoided-crossing cluster with 
Husimi structures that appear as mixtures between the 
initial Husimi structures of the 120th and 121st levels. 
The 121st state, however, leaves with a Husimi struc- 
ture nearly identical to the initial structure of the 119th 
state. These observations are consistent with the slope 
exchanges in Fig. [5] This mixing causes the Husimi dis- 
tributions of the 119th and 120th levels to delocalize after 
the avoided-crossing cluster. Thus, in contrast to sharp 
avoided crossings, broad avoided-crossing clusters play a 
significant role in modifying the aggregate Husimi struc- 
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m=18 m=19 




FIG. 9: [Color online] Husimi structure exchange through the 
sharp avoided crossing shown in Fig. [3b. The left and right 
columns show the Husimi distributions of the 18th and 19th 
levels, respectively. The harmonic oscillator momentum IT is 
on the vertical axis and the oscillator position <T> is on the 
horizontal axis. During the structure exchange, probability 
flows continuously from the integrable region to the chaotic 
region for the 18th eigenstate (and vice- versa for the 19th). 

ture as the coupling strength a increases. In particular, 
such clusters appear to mix the Husimi structure of in- 
dividual eigenstates. 



Localization and a.- dependence of Avoided Crossings 

We have observed that sharp avoided crossings involve 
an interaction either between two chaotic levels or be- 
tween a level localized primarily in an integrable region 
(with only slight localization in the chaotic region) and 
one localized primarily in a chaotic region. However, 
broad avoided crossings such as those in Fig. [5] tend to 
involve interactions between levels with significant local- 
ization in both integrable and chaotic regions. 

We have also seen that broad avoided crossings become 
more prevalent as a increases (see, for example, Fig- Sfc-) ■ 



m=119 m=120 m=121 




C F 



FIG. 10: [Color online] Mixing of Husimi structures through 
the avoided-crossing cluster of Fig. [5] The left, center, and 
right columns show the Husimi distributions of the 119th, 
120th, and 121st levels, respectively. The axes are as in Fig. [9] 



Because clusters of broad avoided crossings lead to mix- 
ing of the Husimi structures of the participating levels, 
one expects that mixing among eigenstates should also 
become more prominent for larger values of a. Addition- 
ally, the a-step size required to resolve avoided crossings 
(i.e., to confirm that they are avoided crossings rather 
than actual level crossings) appears to be correlated with 
the initial localization of the Husimi distributions of the 
participating energy levels. Excluding avoided crossings 
between two chaotic states, the maximum step size that 
is sufficient for resolution seems to increase as the lo- 
calization of the participating levels in similar regions 
of the classical SOS increases. Because adjacent eigen- 
states with energy-level curves that are nearly parallel 
over some range of a can be interpreted as an extremely 
long-range avoided crossing (with a very large step size 
sufficient for resolution), they provide a limiting case of 
the above situation. Such examples arise when the par- 
ticipating eigenstates are localized to almost the same 
extent in the same integrable region. 

We have observed a related phenomenon in the rate 
at which the doubly degenerate eigenvalues at a = 
split (see Fig. [2]). For example, when a lies strictly be- 
tween and the first sharp avoided crossing in Fig. [2 
the 17th and 18th levels are regular (with Husimi struc- 
ture similar to panel A of Fig. [9|), whereas the 19th and 
20th levels are chaotic (with Husimi structure similar to 
panel D of Fig. [5]). The eigenvalues of the two chaotic 



9 



states clearly diverge from each other faster than those 
of the two regular states. This observation is reminis- 
cent of that discussed above for avoided crossings except 
that the eigenenergies are degenerate at some value of 
a rather than nearly degenerate. Indeed, suppose that 
one is starting from the point of closest approach of an 
avoided crossing. As a increases, a broad avoided cross- 
ing between two states initially localized primarily in in- 
tegrable regions will display a much smaller rate of diver- 
gence from near-degeneracy than a sharp avoided cross- 
ing between two chaotic states. We have also observed 
that degenerate eigenstates always split into either two 
chaotic states or two regular ones and that degenerate 
chaotic states always seem to split at a faster rate than 
degenerate regular ones. 



Signatures of Chaos 

Because ergodicity and exponential divergence of 
phase-space trajectories can be used to characterize clas- 
sical chaotic systems, it has been suggested that dereal- 
ization in the Husimi distributions of a quantum system 
is a possible signature of chaos in its classical counter- 
part [491 ]. This has been quantified and studied in nu- 
merical investigations HH and is also germane to 
the system investigated here, as the fraction of phase 
space with chaotic dynamics in the classical system ([T]) 



increases with L 12fl . Consequently, the derealization 



and mixing of Husimi structures, which both become 
more prominent as the prevalence of broad avoided cross- 
ings increases with L, seem to be signatures of the chaos 
in the corresponding classical system. 



CONCLUSIONS 



the a-step size required to resolve an avoided crossing 
is correlated with the extent to which the participating 
states are localized in the chaotic and integrable portions 
of phase space. 

As the size of the non-interaction region increases, 
the avoided crossings broaden and their density in- 
creases. This, in turn, leads to an increase in the num- 
ber of avoided-crossing clusters, in which the partici- 
pating energy-level curves do not fully exchange slopes 
(as a function of a). We showed numerically that such 
avoided-crossing clusters mix the Husimi structures be- 
tween participating states rather than exchange them 
fully as in the sharp avoided crossings. Such mixing 
tends to promote derealization in the eigenstates as the 
coupling strength is increased. This causes a nontriv- 
ial modification in the aggregate Husimi structure as the 
coupling strength is varied, in contrast to the preserva- 
tion of the aggregate structure that is characteristic of 
sharp avoided crossings. Consequently, as the length of 
the non-interaction region increases, one observes an in- 
creasing amount of mixing in the aggregate Husimi struc- 
ture as a function of the coupling strength. This is a 
signature of the dynamics of the corresponding classical 
system, for which the chaotic portion of phase space in- 
creases with the size of the non-interaction region. Thus, 
the appearance of broad avoided crossings, eigenstate de- 
localization, and the mixing of Husimi structures seem to 
be signatures of chaos in the quantum system. Our nu- 
merical computations also suggest that the sharpness of 
avoided crossings is positively correlated with the extent 
to which the participating Husimi structures are local- 
ized in different regions of phase space. Therefore, the 
dynamics of avoided crossings in quantum systems seems 
to be strongly related to the chaotic dynamics of their 
classical counterparts. 



In this paper, we investigated a model system with 
mixed regular and chaotic dynamics that consists of the 
quantization of a one-dimensional free particle on a ring 
coupled to a one-dimensional harmonic oscillator. By ex- 
amining the eigenenergies as a function of the system pa- 
rameters (the coupling strength a and the relative sizes of 
the interaction and non-interaction regions) and comput- 
ing Husimi distributions, we studied the quantum signa- 
tures of the mixed dynamics. We identified key integrable 
and chaotic structures of Husimi distributions by com- 
paring them with the corresponding classical surfaces of 
section. For example, we examined sharp avoided cross- 
ings between states localized in chaotic regions and those 
localized in integrable ones and demonstrated numeri- 
cally the concomitant complete exchange of their Husimi 
structure 2lJ, |5lJ, |53( . We also showed that an avoided 
crossing between two mixed states tends to be broader 
than that between a predominately regular state and a 
predominately chaotic one. Furthermore, we found that 
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Appendix I: The Hamiltonian Matrix 



By the definition of the uncoupled basis, 



Let Hi — Tipart be the Hilbert space for a free par- 
ticle traversing a ring of length 2 + L and H2 = H-osc 
be the Hilbert space for a harmonic oscillator. De- 
fine |n)i = \ip^ art ) and El, respectively, to be the 
nth eigenstate and corresponding nth eigenenergy for 
the particle. We calculate the coordinate-space projec- 
tions {^(3)^-7 }~=i of {|n) l5 1 - n)i}- =1 and their 
eigenenergies from the Schrodinger equation with peri- 
odic boundary conditions: 



dq 



^ n art (q + l(2 + L))=r n art (q),l€Z. 
The (normalized) solutions are 

1 { ( 2imi 



V2+X 



exp 



2 + L 



El = 5- , n e Z - {0} . 

(2 + L) 2 



(0) 



The ground-state energy of the particle is = 
Air 2 /[(2 + L) 2 }. 

Define |fc) 2 = \^J% SC ) an d &k to ^ e tne ^ tn e ig en_ 
state and corresponding fcth eigenenergy of the oscillator. 
Here, Ei = k + 1/2, so that E 2 = 1/2 is the ground-state 
energy 11]. Using the operator definitions ((U), we write 



4?r 2 n 2 / 1 \ 

(Ei)„„, = J^—rf Snn ' ' ( Ea )fcfc' = (^ fc + 2 J Skk ' ■ 

(9) 

The coordinate-space projections for the free-particle 
eigenstates yield 

(Wx) rm , = r L V n art (qT^ rt (q)x(q)dq 
Jo 



r n art (qy€> rt (l)dq- 



Hence, with Eq. ([6|), we obtain 

I o — 7 — ~ — 7\ I —i + ie 2+i " ) if n + n' , 

I 2TZ if n = »' . 

(10) 

Finally, the creation/annihilation operator identities 

a^\k) 2 = VkTT\k + l) 2 , a\k} 2 = Vk\k - 1)2 (11) 
give 

(w 2 ) fefe , = (vFTT4. fe '+i + VF^.fe'-i) . (12) 

Appendix II: The Husimi Distribution 



V2 



($ - in) 



V2 



($ + in) . 



which are, respectively, the creation and annihilation op- 
erators for the harmonic oscillator ll|, |44| . The Hamil- 
tonian (J5J) becomes 

H = ha + - £ - (at + a) *(?). (7) 

The matrix representation of |7} in the uncoupled basis 
\n)i ® I As) 2 [with integer indices n € (— 00,— 1] U [l,oo) 
and k £ [0, 00)] is 



H = E 1 ®I + I®E 2 - aWi 
where I is the identity matrix and 



W 2 



(8) 



( E l)nn' 



-P 



a<a+- 



(E 2 ) fcfc , = 

(W a ) nn , = H*M '<'/ 
(W 2 ) fefe , = / k 



k' 



V2 



(a f + a) 



k' 



The Husimi distribution H^(jj, q, $,n) of a state \tl>) 
of the two-dimensional quantum-mechanical system ([5]) 
is 



^(p,g,$,n) = 7V|(% ? -^ in |)l^)l 2 



(13) 



where IV^pgiri)) is a coherent state localized around 
(p, q, $, II) and iV is a normalization constant. We con- 
struct a coherent state for the coupled system as the 
product 

l%s,g,*,n)) = \^{p,q))i ® IV , (4,fl)>2 , (14) 

where \ip(p,q))i is the coherent state for the free parti- 
cle and |'0($,n))2 is the coherent state for the uncoupled 
harmonic oscillator. The latter is given by [ll[ 



l%>,n)>2 = e 



= p -M* 2 +ri 2 ) 



($ + iH) 



fc=0 



-\kh 



(15) 



Because the uncoupled particle system is (2 + L)-periodic 
in q and p £ K, the phase space is cylindrical. Using the 
procedure of Ref. [46j to define 1^(^,9,4,11)) f° r this topol- 
ogy, we require that the coherent state IV'feg))! satisfies 



(q\4>( P ,q))i = {q + K 2 + L)\^ iP ,q))i, lez. 



(16) 
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One can construct the coherent states \tp(eq))i from 
Euclidean coherent states \ip(p,q))i (which are Gaussian 
wavefunctions) by wrapping them around the cylinder 
and summing overlapping portions. This yields 

oo 

= C* E (9 + K2 + £)fe)>t, (17) 

l — -co 

which satisfies (fTo) and converges because (q + 1(2 + 
L)\?p(p,q)}i is Gaussian. In Eq. (fTT)) . the quantity C is 
a normalization constant to be determined by the condi- 
tion {lp(p.q)\lp(p.q))l = 1. 

The projection of \ip(p.q))i onto the uncoupled particle 
basis {\n)i, | - n)i}£Li is 

r2+L 

HV>(p,9))i = / (n|g)i(g|^(p,g))i d<7 
Jo 

oo ,-2+L 

= ch (n\q + l(2 + L)) 1 (q + l(2 + L)\i, m ) 1 dq. 

i=-oo J ° 

Using Eq. © and the (2 + L)-periodirity of ^l art (q) = 
(q\n)i, we obtain 



To compare Husimi distributions with classical Poincare 
surface of sections, we take p = (2E — $ 2 — II 2 + 2al>) 1 / 2 
(the eigenstate has energy E) and an appropriate value 
of q (with, for example, q — 1 for a section in the center 
of the interaction region). The value for p arises as a slice 
along the classical energy shell, and the value for q cor- 
responds to the choice of Poincare SOS for the classical 
system pj]. 



{n\q)i{q\^{p,q))idq. (18) 



The coordinate-space projection of the Euclidean coher- 
ent state is 



= ( - J exp (-- (q - qf + ip (q - | 



Using the expression for (q\n)i from ([6]), we evaluate ([T8 
and obtain 



H^(p,g)>i = C * exp ^-i(ri-p) 2 -ig(n-| 

(20) 

With the normalization condition (^(p.^lV'Cp,?))! = 1) wc 
determine from (l20t that 



(21) 



[1 

[2 

[3. 

[4] 
[5] 

[6 

[-: 

(19) [9 

[10 
[11 



Thus, if a state \ip) is expressed in the uncoupled ba- 
sis as \ip) = J2n 1 ioJ2'k=o a nk\n)i <8> \k) 2 , we obtain from 
(fHj) . (TT5D, and ® that 



^(p,g,$,IT) 



g(p) 
V2 



n^O k=0 



$ + ti 



x exp 



$ 2 + n 2 + (n - py 



(-§) 



(22) 



In practice, we truncate the sum in Eq. |22J in order 
to compute the Husimi distribution for eigenstates cal- 
culated using a truncated Hamiltonian matrix for (f5]). 
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